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ABSTRACT 

An  ocean  model  response  lo  Hurricane  Ivan  (2004)  over  Ihe  northwest  Caribbean  Sea  and  Gulf  of  Mexico  is 
evaluated  to  guide  strategics  for  improving  performance  during  strong  forcing  events  in  a  region  with  cn 
crgctic  ocean  fcaiurcs  with  ihe  uliimaic  goal  of  improving  coupled  tropical  cyclone  forecasts.  Based  on  prior 
experience,  a  control  experiment  is  performed  using  quasi-optimal  choices  of  initial  ocean  fields,  atmospheric 
forcing  fields,  air  sea  flux  paramelerizalions,  vertical  mixing  parameterization*,  and  both  horizontal  and 
vertical  resolutions.  Alternate  experiments  are  conducted  by  altering  one  single  model  attribute  and  com¬ 
paring  the  results  to  SST  analyses  and  moored  ADCP  current  measurements  to  quantify  the  sensitivity  to  that 
attribute  and  identify  where  to  concenlrate  model  improvement  efforts.  Almospheric  forcing  thai  does  no! 
resolve  the  eye  and  evewall  of  the  storm  (scales  >10  km)  substantially  degrades  the  ocean  response.  Ordering 
other  model  attnhutes  from  greatest  to  least  sensitivity,  ocean  model  initialization  with  regard  to  the  accuracy 
of  upper-ocean  temperature-salinity  profiles  along  with  accurate  location  of  (Kean  currents  and  eddies  is  the 
most  important  factor  for  ensuring  good  ocean  model  performance.  Ocean  dynamics  ranks  second  in  this 
energetic  ocean  region  l>ceause  a  one-dimensional  ocean  model  fails  to  capture  important  physical  processes 
that  affect  SST  cooling  Wind  stress  drag  coefficient  paramelerizalions  that  yield  values  exceeding  2.5  x  10  •' 
at  high  wind  speeds  or  that  remain  <2.0  X  10  y  over  all  wind  speeds  reduce  the  realism  of  wind-driven  current 
profiles  and  have  a  large  impact  on  hoth  SST  cooling  and  the  heat  flux  from  ocean  to  atmosphere  Turbulent 
heat  flux  drag  coefficient  parameterization*  substantially  impact  the  surface  heat  flux  while  having  little 
impact  on  SST  cooling,  which  is  primarily  controlled  by  entrainment  at  the  mixed  layer  base.  Vertical  mixing 
paramctcrizations  have  a  moderate  impact  on  SSI  cooling  hut  a  comparatively  larger  impact  on  surface  heal 
flux.  The  impacts  of  altering  the  horizontal  and  vertical  resolutions  are  small,  with  horizontal  resolution  of 
—  10  km  and  vertical  resolution  of  —  10  m  in  the  mixed  layer  being  adequate.  Optimal  choices  of  all  attributes 
for  simulating  the  ocean  response  to  Ivan  are  identified 


1.  Introduction 

Coupled  ocean-atmosphcrc  general  circulation  models 
have  become  increasingly  important  for  tropical  cyclone 
(TC)  forecast  guidance  at  operational  prediction  centers. 
Implementation  and  advancement  of  coupled  TC  forecast 
models  such  as  the  Geophysical  Fluid  Dynamics  Labo¬ 
ratory  (GFDL)  model  and  more  recently  the  Flurrtcane 
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Weather  Research  and  Forecast  (HWRF)  model  have 
significantly  improved  track  forecasts  by  the  National 
Hurricane  Center  (e.g..  Bender  et  al.  2007).  However, 
little  improvement  has  been  realized  for  intensity  fore¬ 
casts.  Given  that  the  ocean  provides  the  thermal  energy 
for  intensification,  errors  and  biases  in  the  ocean  compo¬ 
nent  of  coupled  TC  forecast  models  may  contribute  to  the 
low  skill  of  intensity  forecasts. 

When  atmospheric  conditions  arc  favorable.  TC  in¬ 
tensification  often  occurs  as  a  storm  passes  over  regions 
with  high  upper-ocean  heat  content.  This  is  particularly 
true  for  potentially  dangerous  rapid  intensification.  The 
impacts  of  the  Loop  Current  (LC)  and  warm-core  anti 
cyclones  in  the  Gulf  of  Mexico  (GOM)  are  documented 
for  Hurricanes  Gilbert  (1988)  and  Opal  (1995)  (Jacob 
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cl  al.  2CXX):  Hong  cl  al.  2000;  Shay  ct  al.  21XX);  Jacoh  and 
Shay  2003),  and  also  for  Hurricanes  Katrina  and  Rita 
(2005)  (Scharroo  el  al.  2005;  Sun  cl  al.  2006.  Shay  2009; 
Jaimes  and  Shay  2009,  2010).  Similar  results  have  been 
obtained  in  the  western  Pacific  for  Typhoon  Macmi  (2003) 
(Lin  et  al.  2005)  and  have  also  been  obtained  in  statistical 
studies  of  multiple  storms  (Wada  and  Usui  2007;  Lin  et  al. 
2008).  Numerical  models  have  reproduced  the  positive 
impacts  of  high  heat  content  on  intensification  (Schade 
and  Emanuel  1999;  Bender  and  Ginis  2000;  Hong  et  al. 
2000;  Emanuel  et  al.  2 004;  Lin  et  al.  2005;  L.  Wu  et  al.  2005; 
C-C  Wu  et  al.  2007).  Low  ocean  heal  content  can  inhibit 
intensification,  possibly  contributing  to  the  weakening  of 
both  Ivan  (Walker  et  al.  2005)  and  Rita  (Sun  et  al.  2006; 
Shay  and  Uhlhom  2008;  Jaimes  and  Shay  2009.  2010)  as 
they  passed  over  cold-core  cyclones  in  the  COM. 

To  correctly  forecast  intensity  evolution,  the  ocean 
component  of  coupled  forecast  models  must  accurately 
predict  the  rale  and  pattern  of  SST  cooling  relative  to  the 
eye  of  the  hurricane.  However,  ocean  models  have  not 
been  thoroughly  evaluated  for  this  purpose.  The  present 
study  evaluates  an  ocean  model  response  to  Hurricane 
Ivan  (2(X)4)  over  the  northwest  Caribbean  Sea  anti  GOM. 
Evaluation  is  difficult  because  ocean  general  circulation 
models  incorporate  a  large  suite  of  numerical  algorithms 
and  subgrid-scale  para meterizat ions  of  processes  not  ex¬ 
plicitly  resolved  by  the  model  grid.  Surface  fluxes  of 
momentum,  heal,  and  mass  (evaporation)  must  also  be 
parameterized.  The  ocean  slate  must  be  accurately  ini¬ 
tialized.  while  the  impacts  of  the  ocean  outside  the  model 
domain  must  be  provided  with  specified  open-ocean 
boundary'  conditions.  Model  performance  may  also  be  af¬ 
fected  by  the  horizontal  and  vertical  resolutions  of  the  grid. 
ITic  relative  impacts  of  these  attributes  on  predicting  the 
SST  response  arc  poorly  known,  and  it  is  not  clear  where 
the  greatest  effort  toward  improving  model  performance 
should  be  targeted.  Model  improvement  will  he  difficult 
hecause  errors  and  hiases  in  ocean  simulations  arise 
simultaneously  from  all  of  these  factors.  For  example, 
modifications  of  the  vertical  mixing  parameterization  that 
improve  the  upper-ocean  response  may  actually  be  com¬ 
pensating  for  errors  and  hiases  resulting  from  surface  llux 
parametcrizations.  A  credible  model  improvement  effort 
will  require  a  thorough  understanding  of  how  model  nu¬ 
merics  and  parametcrizations  impact  the  predicted  upper- 
ocean  response  to  TC-  forcing. 

The  overarching  goal  of  this  analysis  is  to  determine  how 
to  best  invest  our  collective  efforts  toward  improving  (Kean 
model  performance.  Within  this  context,  sensitivity  to  eight 
individual  model  attributes  is  examined:  1)  horizontal  res¬ 
olution,  2)  vertical  resolution,  3)  vertical  mixing  and  vis¬ 
cosity  parametcrizations.  4)  wind  stress  drag  coefficient 
parameterization,  5)  turhulcnt  heal  flux  drag  coefficient 


parameterization,  6)  atmospheric  forcing  resolution,  7) 
(Kean  model  initialization,  and  8)  (Kean  dynamics  (three 
versus  one-dimensional).  Multiple  experiments  are  per¬ 
formed,  with  one  identified  as  the  control  experiment  and 
the  others  identical  to  it  except  for  altering  one  single 
attribute.  Sensitivity  to  each  attribute  is  quantified  and 
their  relative  levels  of  importance  are  ranked  without  the 
complicating  influence  of  atmospheric  feedhacks  that 
would  he  present  when  using  a  coupled  model.  Analysis  is 
performed  on  three  fields  based  on  their  expected  level  of 
importance  regarding  storm  intensity.  In  addition  to  the 
obvious  choice  of  SST,  ocean-atmosphere  turbulent  heat 
flux  (latent  plus  sensible)  averaged  within  specified  radii 
of  the  storm  center  along  with  upper-ocean  velocity  pro 
files  are  considered.  Heal  flux  is  directly  associated  with 
impacts  of  the  ocean  on  intensity  while  accurate  repre 
sentation  of  velcKity  profiles  is  necessary  for  the  model  to 
reproduce  the  shear-driven  turhulence  at  the  ocean 
mixed  layer  (OML)  base  that  is  primarily  responsible  for 
SST  cooling. 

In  addition  to  this  sensitivity  analysis,  we  assess  the 
overall  realism  of  the  simulated  ocean  response  to  Hur 
ricane  Ivan  within  the  limits  of  availahlc  observations, 
specifically  an  SST  analysis  generated  from  satellite  and 
in  situ  observations  and  velocity  profiles  obtained  from 
moored  acoustic  Doppler  current  profiler  (ADCP)  mea¬ 
surements.  Unfortunately,  targeted  aircraft  observations 
that  could  have  provided  simultaneous  subsurface  pro¬ 
files  of  temperature,  salinity,  and  currents  before,  during, 
and  after  Ivan  were  not  available,  thus  limiting  the  extent 
of  the  evaluation  that  could  be  performed 

Section  2  describes  the  ocean  model,  forcing  fields, 
and  observations.  Section  3  summarizes  the  model  exper¬ 
iments  and  analysis  pnKedures.  The  evaluation  and  sensi¬ 
tivity  analysis  are  presented  in  the  next  three  sections,  with 
section  4  focusing  on  the  SST  response  pattern,  section  5 
emphasizing  thermal  fluxes  from  the  (Kean  to  the  atmo 
sphere  during  Ivan,  and  section  6  addressing  the  ocean 
dynamics,  specifically  wind  driven  upwelling  and  forced 
upper-ocean  velocity  fluctuations  that  are  important  for 
shear-driven  OML  entrainment.  ( '(including  remarks  arc 
presented  in  section  7. 

2.  Model  and  observations 

a.  HYCOM 

The  Hybrid  Coordinate  Ocean  Model  (HYCOM)  is  a 
primitive  equation  ocean  model  that  uses  a  hybrid  vertical 
coordinate  designed  to  quasi-optimallv  resolve  vertical 
structure  throughout  the  ocean.  This  coordinate  system  is 
isopycnic  in  the  stratified  ocean  interior,  but  dynamically 
transitions  to  level  coordinates  near  the  surface  to  provide 
resolution  in  the  surface  mixed  layer  and  to  either  level  or 
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terrain- following  (a)  coordinates  in  the  coastal  ocean.  This 
strategy  enables  HYCOM  to  use  advanced  turbulence 
closures  for  vertical  mixing  and  also  to  be  used  as  both  a 
coastal  and  open-ocean  model  while  retaining  the  advan¬ 
tages  of  isopycnic  coordinates  in  the  stratified  ocean  in¬ 
terior.  Model  equations  and  an  initial  evaluation  of  the 
hybrid  vertical  grid  generator  are  presented  in  Bleck 
(2002).  Subsequent  evolutions  and  further  evaluations 
of  the  model  are  summarized  in  Chassignct  et  al.  (2003, 
2007)  and  Halliwell  (2004). 

/>.  Model  initialization 

All  experiments  except  one  arc  initialized  with  fields 
provided  by  a  data-assimilative  hindcast  that  employs  the 
Navy  Coupled  Ocean  Data  Assimilation  (NCODA)  opti¬ 
mum  interpolation  system  (Cummings  2005)  to  assimilate 
satellite  altimetry  and  SST  along  with  in  situ  observations. 
Altimetry  assimilation  ensures  that  major  ocean  features 
such  ns  the  LC  and  associated  warm-core  anticyclones 
and  cold-core  cyclones  arc  initialized  in  their  correct  lo¬ 
cations.  This  is  especially  important  in  the  GOM  and 
Caribbean  Sea  due  to  the  limited  availability  of  in  situ 
observations.  The  NCODA  system  uses  the  Cooper  and 
Haines  (1996)  algorithm  to  provide  a  downward  projec¬ 
tion  of  anomalous  temperature  and  salinity  profiles  as¬ 
sociated  with  anomalous  sea  surface  heights  (SSHs). 
Halltvvell  et  al.  (2008)  evaluated  the  same  GOM  NCODA 
hindcast  product  regarding  ocean  model  initialization  for 
several  storms  including  Ivan,  noting  a  consistent  cold  bias 
that  was  not  as  severe  before  Ivan  as  before  other  storms. 

To  assess  the  sensitivity  to  initialization,  one  expert 
ment  is  initialized  from  a  nonassimilative  GOM  simu¬ 
lation  where  the  ocean  features  are  unconstrained  by 
observations.  Significant  differences  exist  in  the  LC  path 
and  in  the  locations  of  warm-core  anticyclones  and  cold- 
core  cyclones  between  the  two  initialization  products 
(Fig.  I ).  The  two  cold-core  cyclones  present  in  the  data- 
assimilative  initialization  (Fig.  1,  red  arrows)  that  are  crit¬ 
ically  important  for  simulating  the  SST  response  pattern 
forced  by  Ivan  (Halliwell  et  al.  2008)  are  not  present  at 
the  same  locations  in  the  non-assimilativc  initialization. 
For  all  experiments,  the  outer  model  product  used  for  ini¬ 
tialization  also  provides  the  required  open-ocean  boundary 
conditions  during  the  runs. 

c.  Surface  atmospheric  forcing 

Realistic  atmospheric  forcing  fields  that  resolve  the 
inner-core  structure  of  hurricanes  arc  required  to  per¬ 
form  meaningful  ocean  model  evaluation.  Atmospheric  or 
coupled  occan-atmosphcrc  models  run  globally  by  major 
operational  centers  presently  have  horizontal  resolutions 
that  are  far  too  coarse.  The  most  realistic  available  choice 
is  the  Coupled  Ocean-Atmosphere  Mesoscale  Prediction 
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FlG.  t.  SSH  maps  al  0000  UTC  10  Sep  2(XW  from  inilial  model 
fields  provided  by  (lop)  a  dala-assimilalive  GOM  hindcasl  for 
experiments  GOM1-GOMI3  plus  GOM15  and  (bottom)  a  non- 
asstmilaiive  GOM  Simula  lion  for  experimeni  GOM  14.  ITie  red  arrows 
poini  to  the  l  wo  cold -core  cyclones  where  b  an  forced  large  cooling 


System  (CO  AM  PS;  Hodur  1997)  with  fields  available  at 
27-km  resolution.  Since  even  this  resolution  is  UK)  coarse, 
we  use  the  objectively  analyzed  10  m  vector  wind  fields 
(H*W1ND;  Powell  et  al.  1998).  from  the  National  Oceanic 
a  nd  A  t  mosphe  ric  A  dm  inistrat  ion/  Atlantic  Ocea  nogra  phic 
and  Meteorological  I^aboratory/Hurricane  Research  Di¬ 
vision  (NOAA/AOMIVHRD).  to  improve  the  inner-core 
structure.  Briefly,  these  fields  are  prtxluced  by  the  objec¬ 
tive  analysis  of  aircraft  observations  combined  with  other 
available  meteorological  observations.  The  H*W1ND 
fields  accurately  represent  both  the  maximum  wind  speed 
and  the  radius  of  maximum  wind  speed  (/(MAX),  the  latter 
being  approximately  45  km  for  Ivan.  The  low  resolution  of 
the  COAMPS  product  results  in  an  R\\ax  that  is  about 
twice  as  large  as  observed  and  a  maximum  wind  that  is 
reduced  by  40% -50%. 

Ideally,  higher-resolution  atmospheric  forcing  fields 
from  a  fine-resolution  coupled  model  such  as  HWRF 
should  be  used.  However,  no  such  product  presently  exists 
for  Ivan  where  the  storm  both  follows  the  correct  path  and 
maintains  the  correct  intensity  with  sufficient  accuracy 
at  all  times.  This  is  especially  true  for  studying  the  upper 
ocean  velocity  response  for  reasons  described  in  section  2d. 
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FIG.  2.  Map  of  10-m  wind  speed  in  the  eastern  GOM  on  15  Sep 
from  the  blended  COAMPS-H*WtND  forcing  product  Hie  white 
circle  illustrates  the  radius  of  maximum  wind  =  45  km.  l*he 

magenta  (1.4KMAX)  and  black  (3.0KMAX)  circles  illustrate  the  av¬ 
eraging  domains  for  the  air-sca  thermal  heal  flux  sensitivity  anal¬ 
ysis  in  Fig.  9. 

Wind  stress  and  wind  speed  forcing  fields  are  gener¬ 
ated  following  the  procedure  described  in  Halliwell  et  al 
(2008)  that  blends  vector  wind  maps  from  H*WTND  with 
output  from  the  1°  Navy  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS)  atmospheric  model.  Im¬ 
proved  fields  arc  generated  for  the  present  study  because 
the  H*W1ND  fields  arc  blended  with  the  higher-resolution 
CO  AM  PS  atmospheric  model  and  because  the  H*WIND 
fields  are  now  available  at  higher  temporal  resolution 
(3-hourly  versus  only  at  the  time  of  research  aircraft 
flights).  The  3-hourly  COAMPS  and  H*WTND  wind 
component  fields  were  temporally  interpolated  to  30-nun 
intervals  and  then  blended  using  Fqs.  (5)  and  (6)  in 
Haiti  well  et  al.  (2008)  with  r,  -  160  km  and  r2  ~  420  km. 
Thus,  the  vector  wind  fields  are  100%  (()%)  H*WIND  < 
160  km  (>420  km)  from  the  storm  center.  One  reali¬ 
zation  of  the  blended  wind  speed  forcing  is  presented  in 
Fig.  2.  Vector  wind  stress  fields  are  calculated  from  bulk 
formulas  using  a  prescribed  drag  coefficient  Ct). 

Other  required  forcing  fields  (air  temperature,  specific 
humidity,  net  radiative  heal  fluxes,  and  precipitation)  are 
provided  entirely  by  the  COAMPS  model.  These  fields 
are  horizontally  interpolated  onto  the  ocean  model  grid 
points  and  then  temporally  interpolated  to  30-min  in¬ 
tervals.  Surface  turbulent  (latent  plus  sensible)  heat  fluxes 
are  calculated  during  model  experiments  using  bulk  for¬ 
mulas.  After  Ivan  makes  landfall,  the  wind  speed  and 
wind  stress  forcing  revert  to  the  3-hourly  COAMPS  fields. 

<7.  Observations 

Cooling  patterns  are  evaluated  against  daily  SST  fields 
generated  by  the  objective  analysis  of  in  situ  observations 


along  with  Advanced  Very  High  Resolution  Radiometer 
(AVHRR)  and  microwave  satellite  observations  onto  a 
0.25°  global  grid  ( Reynolds  et  al.  2007),  hereinafter  referred 
to  as  the  ‘‘blended’’  SST.  Simulated  upper-ocean  velocity 
profiles  arc  evaluated  against  moored  velocity  measure¬ 
ments.  Hurricane  Ivan  passed  directly  ov  er  an  array  of 
moored  ADCPs  deployed  by  the  Naval  Research  Labo¬ 
ratory  over  the  continental  shelf  and  slope  south  of  Mo¬ 
bile,  Alabama,  as  part  of  the  Slope  to  Shelf  Energetics  and 
Exchange  Dynamics  (SEED)  project  (Carnes  et  al.  2008; 
Teague  et  al.  2006).  Six  shallow  ADCP  moorings  were 
deployed;  three  each  near  the  (X)-  and  90-m  isobaths.  Fight 
other  ADCP  moorings  were  deployed  over  the  conti¬ 
nental  slope,  four  along  the  500-m  isobath  and  four  along 
the  1000-m  isobath  near  De  Soto  Canyon.  Teague  el  al. 
(2007)  used  this  dataset  to  study  fundamental  properties  of 
the  ocean  response  to  Ivan  over  this  shelf-slope  region, 
which  represents  a  challenge  for  ocean  models  (e.g.,  trail 
sition  between  open-ocean  and  coastal  regimes). 

Herein,  we  evaluate  the  simulated  velocity  response 
to  Ivan  at  four  of  the  ADCP  moorings  at  differing  lo¬ 
cations  relative  to  the  eye  with  different  bottom  depths 
(see  inset.  Fig.  5).  M wrings  7  and  9  arc  located  oxer  the 
continental  slope  of  the  northern  GOM  near  500-m  depth 
and  were  located  under  the  western  and  eastern  parts  of 
the  eyewall.  respectively.  Moorings  1  and  3  are  located 
over  the  middle  continental  shelf  near  50- m  depth  and 
were  also  located  under  the  western  and  eastern  parts  of 
the  eyewall.  respectively.  The  forced  velocity  response 
differs  markedly  over  small  differences  relative  to  the  eye 
due  to  differences  in  wind  vector  rotation  on  opposite 
sides  of  the  eye.  ITiis  analysis  therefore  requires  that  the 
storm  in  the  atmospheric  forcing  fields  follow  the  true  path 
with  high  accuracy  (<10-km  lateral  error).  No  available 
product  other  than  UNWIND  provides  this  accuracy  to  the 
best  of  our  knowledge. 

3.  Numerical  experimentation 

Fifteen  lion-assimilative  HYCOM  simulations  were 
conducted  to  assess  model  sensitivity  to  eight  attributes 
(Table  I ).  All  experiments  were  conducted  within  a  do¬ 
main  spanning  the  GOM  and  northwest  Caribbean  Sea. 
henceforth  referred  to  as  the  GOM  domain,  that  is  sit¬ 
uated  on  a  Mercator  grid  with  0.04  x  0.04  cos <p  degree 
resolution  (***4  km),  where  </>  is  latitude.  T  he  coastline 
follows  the  actual  land-sea  boundary  but  a  minimum 
water  depth  of  2  m  is  enforced.  Freshwater  input  from 
12  rivers  is  included. 

All  experiments  are  referenced  to  a  control  experi¬ 
ment  (GOM1 )  that  has  the  attributes  listed  in  the  middle 
column  of  Table  I.  The  NCODA  GOM  hindcast  within 
which  it  is  nested  was  conducted  with  20  vertical  layers 
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Table  L  Summary  of  lhc  model  simulations  The  first  column 
I isls  lhe  model  aiinbuics  while  lhc  second  column  lists  lhe  specific 
attributes  of  lheconirol  experimem  GOM1.  The  ihird  column  lists 
all  of  the  alternate  experiment  along  with  lhe  new  attribute 


Model 
all  ribulc 

Conirol  expt 
(GOM1 ) 

Alternate 

expts 

Horizontal 

0.04°  Mercator 

COM2  0.08°  Mercator 

resolution 

Vertical 

26  layers,  4-8  m 

GOM  3:  21  layers. 

resolution 

in  OML 

7.5-15  min  OML 

Vertical 

KPP 

GOM4  31  layers. 

3-5  m  in  OML 

GOM5:  MY 

mixing 

GOM6:  GISS 

Donelan 

GOM7  Powell 

Tm  .  c  KS 

COARE  3.0 

GO  MS:  Large  and  Pond 
CiOM9:  Large  and 

Pond  (capped) 

GOM10:  Shay  and  Jacob 
GOM  1 1 :  Jarosz  et  al. 
GOM  12:  Kara  cl  al. 

Atmospheric 

algorithm 

27- km  COAMPS  + 

GOM  13.  27-km 

forcing 

H*W1ND 

COAMPS  only 

Outer  model 

NCODA  GOM 

GOM14  tree  GOM 

hindenst 

simulation 

Ocean 

three-dimensional 

GOM  15.  one  dimensional 

dynamic's 

on  the  identical  0.04°  mesh.  The  control  experiment  has  26 

layers,  with 

additional  layers  ; 

idded  to  provide  vertical 

resolution  of  4-8  m  in  the  OML  (Fig.  3.  middle  panel) 

The  six  new 

layers  are  assigned  isopvcnic  target  densities 

lighter  than  any  water  present  within  the  GOM  domain  to 
force  them  to  remain  in  the  constant-thiekness  layers  at 
the  top  of  the  water  column.  Vertical  mixing  is  provided  by 
the  A'  profile  parameterization  (KPP)of  Uirgeet  al.  (1994) 
that  has  been  modified  in  HYCOM  to  include  a  bottom 
boundary  layer  parameterization  (Halliwell  et  al.  2(X)9). 
l  he  vector  wind  stress  forcing  held  is  calculated  prior  to 
the  model  run  using  the  Cp  representation  of  Donelan 
el  al.  (2(X)4).  The  surface  turbulent  heat  Hux  and  evapo¬ 
ration  rate  arc  calculated  during  the  model  run  using  the 
default  parameterization  of  the  latent  and  sensible  heat  flux 
drag  coefficients  Cp\  and  Cp$%  specifically  version  3.0  of 
the  Coupled  Ocean- Atmosphere  Response  Experiment 
(COARE  3.0)  algorithm  (Fairall  et  al.  19%:  Andreas  and 
DeCosmo  2002).  Model  fields  are  initialized  on  (KXX)  UTC 
10  September  2004  and  the  simulation  is  run  for  3 
weeks. 

The  remaining  experiments  (GOM2-GOM15)  each 
differ  from  GOM  1  by  altering  a  single  model  attribute  as 
listed  in  Table  1  to  examine  the  sensitivity  of  the  ocean 
response.  The  alternate  horizontal  resolution  experiment 
GO  M2  is  run  in  the  same  domain,  but  on  a  0.08°  Mercator 
mesh  consisting  of  every  other  point  of  the  high -resolution 


grid.  Two  alternate  vertical  resolutions  are  tested,  one 
coarser  (21  layers,  GOM3)  and  one  finer  (31  layers. 
GOM4)  than  the  control  experiment  (Fig.  3;  Table  1 ). 
requiring  1  and  1 1  additional  layers,  respectively,  to  be 
added  to  the  outer  model  fields.  Two  alternate  vertical 
mixing  parametcrizations  arc  evaluated:  the  Mcllor 
Yamada  (MY)  level  2.5  turbulence  closure  (Mcllor  and 
Yamada  1982)  in  GOM5  and  the  Goddard  Institute  for 
Space  Studies  (GISS)  level  2  turbulence  closure  (Canuto 
et  al.  2(X)I,  2(X)2)  in  GOM6  For  the  wind  stress  drag  co¬ 
efficient,  experiments  GOM7  GOM  1 1  use  the  algorithms 
of:  Powell  et  al.  (2003),  Targe  and  Pond  ( 1981 ),  Uirge  and 
Pond  with  CD  capped  at  a  constant  value  above  30  m  s  *. 
Shav  and  Jacob  (2006).  and  Jarosz  et  al.  (2007).  These 
parametcrizations  are  graphed  in  Fig.  4.  One  alternate 
representation  of  Cri  and  Cp s  that  is  built  into  HYCOM 
is  used  in  GOM  12,  specifically  the  Kara  et  al.  (2002)  al¬ 
gorithm.  Coefficients  from  both  algorithms  arc  graphed  in 
Fig.  4  for  a  representative  choice  of  air  temperature  and 
SST.  Experiment  GOM  13  evaluates  the  sensitivity  to  at¬ 
mospheric  forcing  resolution  by  using  CO  AM  PS  forcing 
alone  and  demonstrates  the  necessiiv  of  using  blended 
H*\VIND  forcing  for  the  control  experiment.  GOM  14 
is  nested  in  the  non -assimilative  outer  model  to  eval¬ 
uate  sensitivity  to  initialization.  Finally,  GOM  15  uses 
the  HYCOM  code  configured  as  one-dimensional  ocean 
models  run  independently  at  each  grid  point  to  demon¬ 
strate  the  impacts  of  three-dimensional  ocean  dynamics. 

4.  SST  cooling 

a  Results  from  the  control  experiment 

The  control  simulation  is  conducted  with  the  same 
model  parametcrizations  as  the  experiment  analyzed  by 
Halliwell  et  al.  (2(X)8)  to  document  the  impacts  of  ocean 
model  initialization.  However,  improved  comparison  be 
tween  model  and  observations  is  expected  in  the  present 
analysis  because  of  1 )  the  use  of  improved  blended  forcing. 
2)  upgrades  to  the  model  code,  and  3)  the  use  of  a  more 
optimal  SST  analysis  for  model  evaluation.  As  in  the 
earlier  study,  poststorm  SST  cooling  in  the  control  ex¬ 
periment  is  largest  within  the  two  cold-core  cyclones  ini 
tially  present  in  the  eastern  GOM  (Figs  1  and  5).  The 
comparison  of  ni(xlcl  results  to  an  SST  analysis  derived 
solely  from  microwave  satellite  measurements  performed 
by  Halliwell  et  al.  (2(X)8)  demonstrated  that  the  model 
caused  a  substantial  overcooling  in  both  cyclones.  How¬ 
ever.  a  visual  comparison  of  the  post  Ivan  AVHRR  im 
ages  analyzed  by  Walker  et  al.  (2(X)5)  to  the  post-Ivan 
microwave  SST  fields  presented  in  Halliwell  ct  al.  (2008) 
suggested  that  the  daily  microwave  SSI  analysis  may  have 
underestimated  the  actual  cooling  that  occurred  within 
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FIG.  3.  Upper-ocean  zonal  cross  sections  of  initial  ((KXK)  IJTC  10  Sep  2004)  ocean  lemperature  fields  provided  h\ 
the  data-assimilative  GOM  hindcast  for  the  three  different  vertical  resolutions  tested,  (top)  low  resolution  used  for 
CiOM3,  (bottom)  high  resolution  used  for  COM4,  and  (middle)  medium  resolution  used  for  all  other  experiments 
Flic  fixed  coordinates  near  the  surface  follow  the  HYCOM  convention  of  increasing  in  thickness  with  depth  until 
a  prescribed  maximum  thickness  is  reached  The  range  of  laver  thicknesses  for  each  experiment  is  listed  in  Table  1 


j 


March  201 1 


H  A  I  !  IWI  M  !  1  At, 


927 


Wiftd  Sir#**  0r*Q  CwHleMwIC* 


Fi(i.  4  (lop)  1  lie  wind  stress  drag  coefficients  ('r>  as  a  funciion 
of  wind  speed  and  (bottom)  t lie  talent  and  sensible  heal  flux 
coefficients  Cri  and  Cns  as  a  function  of  wind  speed  provided 
to  model  experiments.  The  six  wind  stress  drag  coefficients  are 
Donclan  (black);  Powell  (red);  1  arge  and  Pond  (green);  Large  and 
Pond  capped  (blue  above  30  ms1);  Shay  and  Jacob  (magenta), 
which  reverts  to  Large  and  Pond  below  10  ni  s  and  Jaros/  cl  al. 
(black  dashed),  which  re  veils  to  Large  and  Pond  below  20  m  s  1 
and  levels  off  al  1.5  X  If)  3  above  47  m  s  1  The  heal  flux  co- 
efficienis  are  the  COARE  3.0  algorithm  (CEL,  black  solid;  CLS. 
black  dashed)  and  Kara  et  al.  (C\  l.  red  solid;  Chs»  red  dashed),  and 
are  calculated  assuming  an  air  temperature  of  2TC  and  an  SST  of 
2 9'C.  Similar  differences  between  these  two  algorithms  are  re¬ 
alized  for  different  choices  of  air  temperature  and  SS  I 

these  cyclones.  The  daily,  0.25°  Reynolds  el  al.  (2007) 
blended  analysis  of  microwave.  AVHRR.  and  in  situ  SST 
was  selected  for  the  present  study  because  it  produced  a 
post-Ivan  SST  that  was  colder  by  *»2°C  in  the  northern 
cyclone  and  *»1°C  in  ihc  southern  cyclone,  more  in 
agreement  with  the  AVHRR  images.  With  this  change, 
the  control  experiment  overcoolcd  with  respect  to  the  SST 
analysis  by  only  0.6°  in  the  northern  cyclone  (Table  2). 
However,  it  still  excessively  overcoolcd  by  4.3nO  in  the 
southern  cyclone.  Because  the  simulated  cooling  over 
the  remainder  of  the  eastern  COM  and  northwest 
Caribbean  Sea  is  generally  close  to  the  observed  mag¬ 
nitude  (Fig.  5),  large  errors  were  generally  confined  to 
the  southern  cyclone. 


The  impacts  of  ocean  features  on  the  cooling  pattern  are 
explored  further  by  graphing  temperature  as  a  function  of 
time  and  depth  at  the  two  locations  illustrated  in  Fig.  5 
(top):  one  located  at  the  eastern  edge  of  the  detached 
warm  anticyclone  and  the  other  located  near  ihc  center  of 
the  southern  cold-core  cyclone.  Both  points  are  located 
about  45  km.  or  one  radius  of  maximum  wind  speed 
(/?max),  cast  of  the  track  with  estimated  from  the 

blended  wind  speed  forcing  field  (Fig.  2).  This  is  the  ap¬ 
proximate  location  relative  to  the  track  where  the  largest 
near-inertial  currents  are  forced  (e  g.,  Price  1981;  Sanford 
et  al.  2007).  These  profiles  arc  calculated  during  model 
runs  by  inserting  synthetic  moonngs  at  these  two  locations 
with  instruments  deployed  at  5  m  vertical  intervals.  Model 
fields  are  sampled  at  each  synthetic  instrument  using  two- 
dimensional  polynomial  horizontal  interpolation  to  the 
mooring  locations  and  then  performing  linear  vertical  in¬ 
terpolation  to  instrument  depths. 

SST  cooling  is  larger  withm  the  cyclone  (10°C)  com¬ 
pared  to  the  anticyclone  (3°C)  because  of  both  the  thinner 
initial  surface  warm  layer  and  the  stronger  wind-driven 
upwelling  (Fig.  6).  To  highlight  the  contribution  of  up- 
welling.  the  terms  of  the  OML  thickness  balance  are  es¬ 
timated: 


where  is  the  entrainment  velocity  and  w\f  is  the  vertical 
velocity  at  the  depth  of  the  OML  base  r  v/  The  two  terms 
on  the  right  side  are  diagnosed  from  the  model  output,  and 
then  Wf  is  estimated  as  a  residual  The  OML  base  is  di¬ 
agnosed  ns  the  depth  where  the  temperature  is  0.5°C 
colder  than  the  surface  (layer  I  )  temperature.  Vertical 
velocity  is  diagnosed  during  model  runs  at  the  central 
depth  of  each  model  layer  at  each  grid  point  following 
HalLiwcll  cl  al.  (2003).  It  is  then  horizontally  interpolated 
to  the  synthetic  mooring  locations  using  two-dimensional 
polynomial  interpolation  and  then  vertically  interpolated 
to  the  OML  base  using  linear  interpolation. 

These  three  terms  reveal  large  differences  in  the  OML 
depth  balance  between  the  warm  ring  and  cold  cyclone 
(Fig.  6).  Prior  to  the  time  the  eye  passes  the  warm-ring 
location,  w A/  is  small  so  that  w/  ^  —ftzxilfti.  indicating  that 
upwelling  makes  little  contribution  to  the  cooling.  After 
Ivan  passes,  OML  penetration  abruptly  ceases  as 
gradually  increases,  peaking  about  0.2  inertial  periods 
(IPs)  or  about  5  h  after  passage  and  then  remaining  pos 
itive  through  0.4  IP.  During  this  time  interval.  iv/:  is  small 
or  negative  (detrainment)  and  upwelling  still  has  little 
impact  on  cooling.  By  comparison,  upwelling  in  the  cold 
cyclone  makes  a  substantial  contribution  to  cooling  both 
prior  to  and  during  storm  passage,  wilh  wf  peaking  at  the 
time  of  passage  and  remaining  positive  for  about  0.25  IP 
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Fit*;.  5.  Map  of  (top  left)  SSH  on  17  Sep  2004  shortly  after  I van  made  landfall  (track  in  hlnck  line).  illiiMraling  lhe 
local  ions  of  ihc  LC.  lhe  deiachcd  warm  ring,  and  the  two  cyclonic  eddies  (red  arrows)  near  the  time  of  maximum 
cooling  in  lhe  eastern  GOM  A  while  rectangle  marks  lhe  location  of  a  SEED  mooring  (Teague  et  ul.  2007)  and  (lop 
right)  the  inset  at  light  illustrates  mooring  locations  with  the  approximate  track  of  l\an  shown  b\  Ihc  red  line  The 
other  two  unmarked  while  dots  in  the  SSH  map  represent  lhe  locations  sampled  by  synthetic  instruments  and  de 
scribed  in  Fig.  6.  The  remaining  panels  present  SST  maps  for  (left)  10  and  (nghl)  17  Sep  2(KV4  for  (lop)  the  Reynolds 
et  al  (2007)  blended  analysis  of  in  situ  observations  plus  infrared  and  microwave  satellite  data  and  (bottom)  from  the 
control  experiment  GOMl. 


thereafter.  The  rate  of  OML  deepening  is  small  for  about 
0.2  IP  prior  to  passage  so  that  ivf.  **  tv V/.  During  this  time 
interval,  upwelling  increases  cooling  by  lifting  colder  wa¬ 
ter  toward  the  OML  base  and  by  maintaining  a  thinner 
OML,  with  the  latter  impact  resulting  in  increased  OML 
entrainment  by  both  wind-  nnd  buoyancy-lorced  TKEs  at 
lhe  surface  and  by  shear  instability  of  the  wind-driven 
horizontal  flow  at  the  OML  base.  A  thin  OML  also  has 
relative  small  thermal  inertia  and  will  cool  more  rapidly 
for  a  given  entrainment  heat  flux.  Following  hurricane 
passage,  entrainment  rapidly  weakens  and  the  upwelling 


primarily  acts  to  raise  the  OML  base.  Beginning  around 
0.2  IP  after  passage  when  \\  \f  is  small,  the  weak  stratifi¬ 
cation  remaining  at  the  OML  base  permits  rapid  OML 
deepening  due  to  large  entrainment  driven  by  the  shear  of 
the  near-inertial  current  fluctuations.  Hntniinment  cooling 
remains  significant  until  0.5  IP  after  passage  and  maximum 
cooling  is  realized  about  0.4  IP  (10  h)  alter  Ivan’s  passage. 

Further  insights  on  the  impacts  Of  upwelling  and  other 
aspects  of  three-dimensional  ocean  dynamics  are  ob¬ 
tained  by  comparing  the  evolution  of  the  temperature 
and  vertical  velocity  profiles  from  lhe  one-dimensional 
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Tahlh  2.  Post -I  van  SST  (  C)  on  17  Sep  2004  wiiliin  ihc  northern 
and  southern  cold-core  cyclones  where  the  largest  cooling  occurred 
from  both  the  observations  and  (he  15  experiments  along  with  SST 
differences  between  each  experiment  and  the  observations.  SST 
values  were  calculated  by  averaging  over  1°  x  1°  boxes  centered  at 
28.4TM,  87.4°W  (northern  cyclone)  and  25.0‘N.  86.8AV  (southern 
cyclone). 


Northern  cold-corc  cyclone  Southern  cold-core  cyclone 


Source 

SST 

(T) 

Diff 

from 

GOMI 

Diff 

from 

obs 

SST 

rc> 

Diff 

from 

GOMI 

Diff 

from 

obs 

Obs 

22.3 

23.9 

GOM  1 

21.7 

0.6 

19.6 

4.3 

COM2 

21.6 

0.1 

0.7 

19.5 

-0.1 

44 

COM3 

21.8 

+  0.1 

-0.5 

19.3 

-0.3 

-4.6 

COM4 

21.6 

-0.1 

-0.7 

20.1 

+  0  5 

-3,8 

GOM5 

22.3 

+  0.6 

0.0 

21.4 

+  1.8 

-2.5 

GOM6 

22.5 

+  0.8 

-0.2 

20.2 

+  0.6 

-3.7 

GOM  7 

22.6 

+0.9 

+0.3 

20.5 

+  0.9 

3.4 

GOM8 

20.7 

-1.0 

-1.6 

19.0 

-0.6 

-4.9 

GOM9 

21.5 

0.2 

-0.8 

19.2 

-0.4 

-4.7 

GOM  tO 

20.2 

—  t  .5 

-2.1 

18.1 

-1.5 

-5.8 

GOM  11 

22.3 

+  0.6 

0.0 

19.8 

+  0.2 

-4.1 

GOM  12 

21.7 

0.0 

0.6 

19.7 

+  0.1 

4.2 

GOM  13 

23.7 

+  2.0 

+  L4 

22.0 

+  2.4 

1.9 

GOM  14 

23.4 

+  1.7 

+  t.l 

26.5 

+  6.9 

+  2.6 

GOM  15 

22.7 

+  1.0 

+  0.4 

22.8 

+  3.2 

LI 

experiment  GOM15  to  those  from  the  control  expert 
ment  at  the  same  two  locations  (Fig.  6).  In  the  absence  of 
wind-driven  vertical  velocity  in  both  the  anticyclone  and 
cyclone,  wf  is  large  prior  to  and  during  passage,  but  then 
decreases  rapidly  wilhin  0.1  IP  after  passage  as  wind 
forcing  relaxes  and  the  OML  approaches  its  maximum 
thickness.  SST  cooling  produced  by  GOMI5  is  smaller 
compared  to  the  control,  2C  versus  3°C  in  the  anticyclone 
and  5°  versus  1()°C  in  the  cyclone  (Fig.  6).  As  discussed 
earlier,  wind-driven  vertical  velocity  is  a  major  contributor 
to  increased  cooling  in  the  cyclone.  The  story  is  different  in 
the  anticyclone  because  the  warm  layer  that  is  present  0.5 
IP  before  the  storm  in  GOM15  is  thicker  than  in  the  con¬ 
trol  experiment,  a  situation  that  likely  contributed  sub 
stantially  to  the  reduced  cooling.  The  different  warm  layer 
thicknesses  occur  because  the  measurement  point  is  lo¬ 
cated  at  the  eastern  edge  of  the  anticyclone  and  the  warm 
layer  there  becomes  thinner  over  the  first  few  days  of  the 
simulation  as  this  feature  propagates  westward.  This  anal¬ 
ysis  demonstrates  the  critical  importance  of  ocean  dynamics 
not  only  because  wind-driven  upwelling  can  contribute 
significantly  to  SST  cooling,  but  also  because  the  upwell¬ 
ing  pattern  can  be  highly  distorted  by  the  background 
vorticity  field  associated  with  preexisting  ocean  features 
(Jaimes  and  Shay  2(X)9,  2010). 

Excess  cooling  produced  by  Ihc  control  experiment  in 
the  southern  cyclone  is  too  large  to  ascribe  to  errors  and 


biases  in  the  blended  SST  product.  To  determine  if  the 
model  initialization  could  have  been  a  factor,  the  ocean 
heat  content  (OHC:  Leipper  and  Volgenau  1972)  refer 
enced  to  the  26°C  isotherm  depth  ( D2t >)  derived  from 
satellite  altimetry,  SST  measurements,  and  ocean  clima¬ 
tology  (Mainelli  ct  al.  2(X)8)  is  compared  to  values  ealeu 
latcd  from  initial  model  fields.  Values  for  the  southern 
cy  clone  arc  estimated  over  a  1°  X  1°  box  centered  at  25.0°N. 
86.8°W.  The  derived  OHC  (D2J  value  is  40  kJ  cm  2 
(48  m)  while  the  initial  model  values  are  25  kJ  cm  2 
(32  m).  This  initial  cold  bias  may  therefore  have  con¬ 
tributed  to  the  overcooling.  However,  the  detailed  hori¬ 
zontal  and  \ertical  structures  of  ocean  features  provided 
by  the  initialization  may  also  have  influenced  the  up 
welling  pattern  and  its  impact  on  cooling.  Although  it  is 
possible  that  the  vertical  mixing  parameterization  could 
contribute  to  the  overcooling,  separating  this  contribu¬ 
tion  from  others  will  be  difficult.  Future  efforts  to  improve 
ocean  model  performance  will  depend  heavily  on  the 
availability  on  detailed,  high  qualit\.  three-dimensional 
ocean  observations  acquired  prior  to,  during,  and  sub¬ 
sequent  to  individual  storms  to  identify  and  separate  the 
contributions  of  these  different  factors. 

b.  Sensitivity  to  model  attributes 

The  sensitivity  of  the  SST  cooling  pattern  is  evaluated 
by  comparing  SST  changes  from  experiments  GOM2- 
GOM15  to  those  from  the  control  experiment  GOMI. 
SST  change  is  calculated  for  the  control  experiment  as 
AT(  -  (73  -  T[)c  and  for  each  alternate  experiment  as 
&Ta  -  (73  -  T])a%  where  T\  (73)  is  the  SST  before  (after) 
Ivan.  The  difference  in  SST  change  between  each  alter¬ 
nate  experiment  and  the  control  experiment  given  by 

AT  -  A T A  A/’{  (2) 

is  analyzed  in  Fig.  7,  where  the  tabulated  RMS  amplitude 
of  A T  quantifies  the  similarity  of  each  alternate  experi 
ment  to  the  control  experiment.  Inset  maps  of  A /are  also 
presented  in  Fig.  7  for  four  of  the  alternate  experiments 
All  statistics  are  calculated  within  the  region  outlined  hv 
the  box  in  these  four  maps,  which  covers  the  eastern 
GOM  and  extreme  northwest  Caribbean  Sea. 

The  simulated  SST  response  is  most  sensitive  to  the 
alternate  initialization  (GOM  14)  with  an  RMS  AT  of 
1.57°C.  This  large  amplitude  is  evident  in  the  inset  map 
(Fig.  7).  l^irge  positive  AT  exists  at  the  locations  of  the 
two  cold-corc  cyclones  because  these  features  were  not 
present  in  the  alternate  initialization.  Substantial  sensi¬ 
tivity  to  initialization  extends  beyond  these  ocean  fea¬ 
tures  and  also  beyond  the  region  directly  forced  by  the 
storm.  Ocean  model  sensitivity  to  initialization  is  there¬ 
fore  likely  to  remain  significant  over  (it  her  ocean  regions 


930 


MONTHLY  WEATHER  R  E  V  I E  W 


VoilJMI  139 


GOM1  (control) 


GOM15  (one-dimensional) 


T  (C),  Warm  Ring  T  (C).  Warm  Ring 


-0.5  0  0.5  1  -0.5  0  0.5  1 

OML  Thickness  Eq.  Warm  Ring  OML  Thickness  Eq.,  Warm  Ring 


OML  Thickness  Eq..  Southern  Cyclone 


Ineriiai  Periods 


OML  Thickness  Eq ..  Southern  Cyclone 


Fl(i.  h.  (top  row)  l  ime  series  of  temperature  over  the  upper  400  m  and  {second  row )  of  the  three  terms  of  I  q.  ( I )  at 
a  location  in  the  detached  warm  ring,  (bottom  two  rows)  The  same  two  plots  at  a  location  in  the  southern  cold  core 
cyclone,  (left)  The  control  experiment  GOM  1  and  ( right)  the  one-dimensional  experiment  GOM 1 5  I  hese  locations 
arc  illustrated  by  white  dots  in  the  top  panel  of  Fig.  5.  The  thick  black  line  denotes  the  OML  base  diagnosed  as  the 
depth  where  the  temperature  reaches  0.VC  below  SST.  T lie  magenta  line  follows  the  20°C  isotherm  in  the  depth  lime 
temperature  plots.  For  the  terms  in  Eq.  (1)  wr  is  the  black  line.  *  u  is  the  blue  line,  and  is  the  red  line. 


with  less  energetic  features.  Even  in  the  absence  of  en¬ 
ergetic  ocean  features,  it  is  still  necessary  to  accurately 
initialize  the  temperature  and  salinity  structure.  The  next 
largest  sensitivity  is  observed  for  the  ocean  dynamics 
(GOM  15)  with  an  RMS  ST  of  1.3 1°C  At  the  other  ex¬ 
treme,  the  SST  response  is  leasl  sensitive  to  vertical  reso¬ 
lution  (RMS  ST  of  0.29°C  for  low  resolution  and  0.26°C 
for  high  resolution),  and  also  to  turbulent  heat  flux  drag 
coefficients  (RMS  S T of  Q.I8°C).  In  the  latter  case,  the  lack 
of  sensitivity  exists  despite  differences  of  about  30%-40% 
in  values  of  Cel  and  Crs  (big.  4).  Although  surface  heat 
flux  is  expected  to  be  sensitive  to  these  different  coefficient 
values  (section  5),  it  has  a  relatively  small  tnfluence  on  SST 
cooling,  which  mostly  arises  from  entrainment  across  the 
OML  base. 


Intermediate  ocean  model  sensitivity  exists  for  the 
remaining  attributes  of  horizontal  resolution,  vertical 
mixing,  wind  stress  drag  coefficient,  and  atmospheric 
forcing  resolution.  The  exception  to  Ihis  is  the  small  RMS 
ST  value  of  ().19°C  realized  for  the  Large  and  Pond  C*> 
capped  at  high  wind  speed,  which  produced  C/>  values  quite 
similar  to  those  of  Donclan  used  in  the  control  experiment 
(Fig.  4).  The  inset  ST  map  for  GOM6  (GISS  vertical  mix¬ 
ing)  along  with  the  mean  S T difference  of  0.2I°C  (Fig.  7) 
demonstrates  that  this  mixing  scheme  produces  less  SST 
cooling  than  KPP  The  same  is  true  for  MY  mixing  (not 
shown).  Paramcteri/ations  of  Cp  that  produce  larger 
(smaller)  values  result  in  greater  (less)  SST  cooling.  The 
inset  map  for  GOM8  (Large  and  Pond  C/d  reveals  the 
enhanced  SST  cooling  produced  by  this  experiment  (Fig.  7). 
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Alternate 

Experiment 

Model  Attribute  Changed 
from  Control  Experiment 

\T  (C) 
Mean 
Diff. 

AT  (C) 
RMS 

GOM2 

Low  horizontal  resolution 

+0.01 

0.47 

GOM3 

Low  vertical  resolution 

+0.05 

0.29 

GOM4 

High  vertical  resolution 

-0.02 

0.26 

GOM5 

MY  mixing 

+0.06 

0.49 

GOM6 

G1SS  mixing 

+0.21 

0.37 

GOM7 

Powell  CD 

+0.13 

0.33 

GOM8 

Large&Pond  CD 

-0.19 

0.48 

GOM9 

Large&Pond  (capped)  CD 

-0.08 

0.19 

GOM10 

Shay&Jacob  CD 

-0.39 

0.75 

GOM11 

Jarosz  of  al.  CD 

-0.04 

0.28 

GOM12 

Kara  ct  al.  Cp  Cs 

-0.13 

0.18 

GOM13 

No  H*WIND 

+0.25 

0.70 

GOM14 

Outer  model 

+0.43 

1.57 

GOM15 

One-dimensional 

-0.17 

1.31 

/ 


\ 


FIG.  7.  Sensitivity  analysis  of  ihc  SST  change  forced  by  Hurricane  Ivan  (track  is  black  line),  as  summarized  by  differences  in  ihc  SSI 
change  (17  -  11  Sep)  calculated  using  Faj.  (2)  between  the  control  experiment  GOM1  and  each  of  rhe  remaining  experiments  GOM2 
GOM  15.  rhe  RMS  amplitude  of  A 7  (°C).  which  represents  lhe  RMS  difference  between  the  SST  changes  forced  b\  the  two  experiments, 
is  tabulated  in  the  third  column.  Four  A  7*  maps  arc  shown  as  examples.  All  RMS  A T values  are  calculated  within  the  rectangular  boxes 
outlined  by  black  lines  in  the  maps  (extreme  NW  Caribbean  Sea  and  eastern  GOM). 


Enhanced  cooling  is  also  produced  by  CiOM. 10  (Shay  and 
Jacob  C n)  while  the  opposite  is  true  for  GOM7  (Powell 
C/>).  This  relationship  results  from  differences  in  the 
strength  of  wind-driven  turbulence'  resulting  from  changes 
in  surface  friction  velocity  and.  also,  from  changes  in  the 
magnitude  ol  wind-driven  flow  that  impact  the  magnitude 
of  both  upwclltng  and  shear  at  the  OML  base. 

A  Taylor  (2001)  diagram  (Fig.  8,  top)  is  used  to  compare 
the  ATa  from  each  alternate  experiment  to  the  A T(  over 
the  subdomain  covering  the  eastern  GOM  and  extreme 
northwest  Canhbean  Sea,  specifically  lhe  box  outlined  in 
the  Fig.  7  inset  maps.  To  construct  this  diagram,  the  mean 
is  first  removed  from  each  A TA  and  \  T(  map,  and  then  all 
maps  arc  normalized  by  the  variance  of  A  T( .  This  diagram 
is  effective  because  each  point  simultaneously  represents 
three  different  metrics.  The  large  black  square  in  the  dia¬ 
gram  is  the  point  that  signifies  a  perfect  comparison 


(correlation  coefficient  of  1.0.  identical  normalized  RMS 
amplitude  of  i.0,  and  RMS  difference  of  zero  between 
maps).  One  plotted  point  quantifies  the  similarity  between 
A I A  and  A7\  .  with  the  RMS  amplitude  of  &7'A  given  bv 
the  radial  distance  from  the  plot  origin  at  the  bottom-left 
corner,  the  RMS  difference  between  the  two  fields  being 
proportional  to  the  linear  distance  between  the  plotted 
point  and  the  large  black  square,  and  the  correlation  co¬ 
efficient  being  a  function  of  the  direction  angle  of  the 
vector  connecting  the  plot  origin  with  the  plotted  point. 

The  greatest  sensitivity  again  results  from  the  alter 
natc  initialization  (GOM  14).  which  has  a  (correlation, 
normalized  RMS  amplitude,  normalized  RMS  differ¬ 
ence)  of  (0.59.  0.58,  0.81)  with  respect  to  the  control 
experiment  (Fig.  8).  The  small  RMS  amplitude  of  0.5S 
is  due  largely  to  the  absence  of  the  two  cyclones  di¬ 
rectly  hit  by  Ivan  and  the  associated  large  cooling.  The 
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one-dimensional  experiment  (COM  15)  has  ii  (correla¬ 
tion,  normalized  RMS  amplitude,  normalized  RMS  dif¬ 
ference)  of  (0,72,  0.97,  0.72)  with  the  control  experiment. 
All  other  alternate  experiments  are  highly  correlated  with 
the  control  experiment  (>0.95,  Fig.  8),  demonstrating  that 
correct  model  initialization  and  three-dimensional  ocean 
dynamics  arc  necessary  to  accurately  predict  the  SST  re¬ 
sponse.  At  the  other  extreme,  the  three  experiments  dis¬ 
playing  the  least  SST  sensitivity  are  GOM12  (alternate 
turbulent  heat  flux  parameterization),  GOM9  (Large  and 
Pond  Co  capped  at  high  wind  speed),  and  GOM15  (Jarosz 
ct  aJ.),  The  dependence  of  RMS  amplitude  on  C/>  is  evi¬ 
dent  as  is  the  several-percent  reduction  in  RMS  amplitude 
produced  by  the  two  alternate  vertical  mixing  choices. 
RMS  amplitude  produced  by  GOMI3  (no  H*WIND 
blending)  is  reduced  by  **25%  due  to  the  highly  smoothed 
representation  of  inner-core  storm  structure  and  the  re¬ 
sulting  reduction  in  maximum  wind  speed. 

Since  these  comparisons  are  substantially  influenced 
by  the  large  cooling  that  occurred  within  the  two  oceanic 
cyclones,  the  focus  now  shifts  to  the  cooling  that  occurred 
in  these  features  (Table  2).  The  largest  underestimate  of 
cooling  relative  to  the  control  experiment  is  produced  by 
GOM14  within  the  southern  cyclone  where  the  post- 1  van 
SST  is  6.9°C  warmer  than  in  GOM1  because  of  the  ab¬ 
sence  of  this  feature.  The  inadequate  resolution  of  storm 
structure  by  the  alternate  forcing  (GOM13)  and  the  ab¬ 
sence  of  three-dimensional  ocean  dynamics  (GOM15) 
each  led  to  reduced  cooling  of  I  .()°-3.2°C  in  both  cyclones 
(Tabic  2).  The  two  alternate  vertical  mixing  choices  re¬ 
duce  the  SST  cooling  relative  to  the  control  by  0.6°-L8°C 
while  Cp  parameteri/ations  that  produce  larger  (smaller) 
values  than  the  Donelan  et  al.  (2004)  choice  over  in¬ 
termediate  to  high  wind  speeds  increase  (decrease)  the 
cooling  by  up  lo  0.9°C  (1»5°C).  Little  difference  in  cooling 
results  from  altering  the  horizontal  and  vertical  resolutions 
or  from  altering  the  turbulent  heat  flux  drag  coefficients 

c.  Evaluation  of  the  SST  response  against 

the  blended  SST  fields 

lo  evaluate  the  realism  of  the  SST  cooling  pattern 
from  all  experiments  against  the  observations,  a  Taylor 


TABt  i  3.  Comparison  of  A  /  maps  (17  10  Sep)  between  the  1 5 

experiments  and  (he  blended  SSI  mean  difference  (column  1) 
RMS  difference  (column  2).  and  Murphy  '-kill  score  S  from  Eq.  (3) 
(column  3).  Skill  scores  are  also  presented  for  u  and  v  profiles  ui 
SEED  mooring  9. 


Ex  pi 

A  /  mean 
diff  (nC) 

A/  RMS 
diff  (rC) 

A/ skill 

score 

t) 

.skill 

score 

0 

skill 

score 

GOMI 

0.24 

1 .33 

0.20 

0.71 

0.77 

COM2 

0.25 

1.29 

0.27 

0.64 

0.73 

COM3 

0.30 

1  40 

0.12 

0.65 

0.7S 

COM4 

(1.22 

1  40 

0.17 

0.63 

0.78 

GOM5 

0.30 

1.26 

0.28 

0.62 

0.61 

GOM6 

0.47 

1.33 

0.07 

0.71 

0.68 

GOM7 

0.3S 

1.33 

0  14 

0  58 

0.60 

GOM8 

0.04 

1.57 

-0.02 

0.47 

0.54 

GOM9 

0.14 

1.47 

0.12 

0.61 

(1.74 

GOMtO 

-0.19 

1.71 

-0.19 

0.45 

0.26 

GOMI  1 

0.06 

1.37 

0.27 

0.66 

0.74 

COM  12 

0.11 

1.34 

0.22 

0.70 

0.76 

GOM 1 3 

051 

1.24 

0.10 

0.42 

0.22 

GOM14 

0.84 

1.51 

0.61 

0.21 

0.65 

GOM15 

0.06 

1.55 

0.04 

-1.10 

0.97 

diagram  analysis  compares  A  T  \ 

and  A  Tc 

maps 

from  all 

15  experiments  to  AT^  maps  calculated  from  the  blended 
SST  analysis  over  the  subdomain  outlined  in  the  Fig.  7 
inset  maps.  In  this  case,  all  maps  arc  first  normalized  by 
the  variance  of  the  A  7'*  map  This  evaluation  is  supple 
mented  by  calculating  both  the  mean  differences  and  the 
RMS  differences  over  the  same  subdomain  (Table  3). 
The  mean  differences  provide  information  that  is  not 
contained  in  the  Taylor  diagram  while  the  RMS  differ¬ 
ences  document  the  dimensional  magnitude  of  the  non- 
dimensional  RMS  differences  that  appear  in  the  diagram. 
The  evaluation  is  also  supplemented  by  the  Murphy 
( 1988)  skill  score  A,  given  by 


r-(Zy) 

■> 

-  W.v /. 

[\  cry  )\ 

w  here  X  and  Y  are  the  mean  v  alues  of  holds  /V  and  V,  <r  v 
and  rrv  are  their  standard  deviations,  and  r  is  the  corre 
lation  coefficient  between  them.  If  the  two  maps  are 
identical,  then  5  -  r\  otherwise.  S  <  r.  This  score  is 
reduced  by  three  factors:  1 )  correlation,  2)  differences  in 


Fic.i  X  (top)  Taylor  diagram  comparing  A7  calculated  using  Eq.  (2)  from  the  control  experiment  to  A  /  from  all  remaining  experiments 
Symbols  representing  the  different  experiments  are  labeled  in  the  legend.  Different  symbol  colors  in  the  diagram  and  in  the  legend 
categorize  the  individual  model  attributes  that  are  varied  horizontal  and  vertical  resolutions  (red)  vertical  mixing  choice  (green),  wind 
stress  drag  coefficient  (blue),  turbulent  heat  flux  coefficient  (cyan),  atmospheric  forcing  resolution  (nugenia)  and  outer  model  choice 
(hlack  diamond)  l'hc  large  black  square  al  a  correlation  of  1  0  and  an  RMS  ampliiudc  of  1.0  represents  a  perfect  comparison,  (bottom) 
I  aylor  diagram  comparing  A  7  calculated  for  each  experiment  to  the  observed  A  7  estimaicd  from  the  Reynolds  SS  I  analysis.  As  shown  in 
the  legend  at  the  bottom,  symbols  and  colors  follow  the  same  conventions  except  that  the  control  experiment  is  included  as  a  hlack  circle 
The  large  hlack  square  again  represents  a  perfect  comparison 
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the  mean,  and  3)  differences  in  the  RMS  amplitude. 
Values  exceeding  0  represent  statistically  significant  skill. 

In  the  Taylor  diagram  (Fig.  8,  bottom),  the  control  ex¬ 
periment  has  a  (correlation,  normalized  RMS  amplitude, 
normalized  RMS  difference)  of  (0.71,  1.20.  0.86)  (sec  Ta¬ 
bic  3  for  the  corresponding  dimensional  RMS  difference). 
The  resulting  skill  score  S  —  0.20  is  significant  but  not  large 
due  to  the  20%  amplitude  overestimate  and  also  to  the 
RMS  difference  \  alues  that  are  nearly  as  large  as  the  RMS 
amplitudes  of  the  two  fields  (Table  3).  Eleven  of  the  15 
experiments  produce  correlation  coefficients  between  0,67 
and  0.72,  with  alternate  model  attributes  primarily  influ¬ 
encing  normalized  RMS  amplitudes.  The  dependence  of 
RMS  amplitude  on  the  value  of  CD  is  particularly  evident 
(Fig.  8,  Table  3).  The  two  experiments  using  the  largest 
values  of  Cf)  (GOM8,  GOM10)  produced  negative  S.  Of 
the  remaining  three  experiments,  the  alternate  initializa¬ 
tion  experiment  GOM14  is  the  largest  outlier,  with  (cor 
relation,  normalized  RMS  amplitude,  normalized  RMS 
difference)  of  (0.41,  0.66,  0.97)  and  insignificant  S.  The 
alternate  forcing  experiment  GOM13  produces  a  sub¬ 
stantially  smaller  normalized  RMS  amplitude  (0.88)  than 
most  other  experiments  along  with  a  slightly  reduced 
correlation  of  0.65  and  small  but  significant  S  =  0. 10.  One- 
dimensional  dynamics  (GOM15)  results  primarily  in  a 
modest  reduction  in  correlation  to  0.61  but  with  a  very 
small  S  =  0.04,  Overall,  the  experiments  overestimate  the 
amplitude  of  the  SST  change  pattern  forced  by  Ivan  by 
^20%  compared  to  the  blended  SST  analysis,  which 
probably  results  to  a  large  extent  from  the  large  over- 
ccxili ng  in  the  southern  cyclone. 

The  mean  differences  listed  in  Table  3  demonstrate  that 
in  the  subdomain  over  which  these  analyses  were  per¬ 
formed,  the  ocean  model  did  not  cool  as  much  as  was 
shown  in  the  blended  SST  fields.  Although  the  model 
tended  to  substantially  overcool  in  the  cyclones,  this  was 
more  than  compensated  for  by  undercooling  over  the  re¬ 
mainder  of  the  domain.  Mean  differences  display  the  larg¬ 
est  sensitivity  to  ocean  model  initialization,  surface  forcing 
resolution,  and  wind  stress  drag  coefficient,  with  more 
(less)  cooling  occurring  for  larger  (smaller)  values  of  Cn. 
Smaller  sensitiv  ity  is  e\  ident  for  the  vertical  mixing  choice, 
with  the  MY  and  G1SS  schemes  producing  slightly  more 
cooling  than  KPP.  Little  sensitivity  is  evident  to  the  vertical 
and  horizontal  resolutions  These  conclusions  also  gener¬ 
ally  hold  for  SST  cooling  within  the  two  cyclones  (Table  2). 

5.  Thermal  energy  provided  by  the  ocean 

to  the  storm 

a.  Estimation  procedure 

Although  the  accuracy  of  SST  forecasts  is  important, 
it  is  the  influence  of  SSI'  cooling  on  surface  turbulent 


heat  (lux  that  most  directly  impactsTC  intensity.  It  is  not 
possible  to  evaluate  SST  cooling  relative  to  Ivan’s  center 
against  observations  as  it  propagates  because  satellite 
sampling  is  restricted  by  the  cloud  cover  and  heavy  pre¬ 
cipitation.  For  these  reasons,  we  do  not  attempt  to  eval¬ 
uate  the  accuracy  of  SST  cooling  following  the  storm. 
Instead,  we  document  the  sensitivity  of  the  surface  tur 
bulcnt  heat  flux  to  model  attributes  b)  averaging  the  flux 
within  a  specified  radius  of  the  storm  center  and  graphing 
it  along  the  storm  path.  The  choice  of  averaging  radius  is 
important  because  it  should  include  the  area  of  the  ocean 
surface  that  is  actually  providing  the  large  majority  of 
thermal  energy  to  the  storm.  Since  this  choice  is  uncertain 
(c.g.,  Shen  et  al.  2002;  Cione  and  Uhlhom  2003),  sensitivity 
analysis  is  performed  for  two  different  radii:  1  4 /?MAX  and 
3/?max>  where  /?Max  's  the  radius  of  maximum  wind. 
These  choices  are  illustrated  in  the  blended  wind  speed 
forcing  map  (Fig.  2),  where  1.4 R\\a\  encloses  the  region 
influenced  by  the  eye  wall,  while  3/?MA\  encloses  most  of 
the  region  influenced  by  hurricane  force  winds.  The  surface 
turbulent  heat  flux  Qj  averaged  over  these  two  radii  is 
graphed  as  a  function  of  storm  latitude  for  all  experiments 
in  Fig.  9.  Although  this  type  of  sensitivity  analysis  is  most 
appropriately  performed  with  a  coupled  model  that  pro 
duces  atmospheric  feedback,  it  is  still  important  to  un 
derstand  how  ocean  model  attributes  alone  influence  the 
oceanic  contribution  to  surface  turbulent  heat  flux  through 
differences  in  SST  cooling 

b.  Sensitivity  analysis 

1‘ he  progression  of  mean  alone  the  storm  path  in  the 
control  experiment  is  strongly  influenced  by  the  oceano¬ 
graphic  features.  Over  the  northwest  Caribbean  and  the 
LC  (south  of  23°N),mean  Q  f  w  ithin  L4/?N1A\  of  the  storm 
center  averages  about  -6(X)  W  m  \  When  the  storm  en¬ 
counters  the  two  cold  cyclones,  SST  cooling  inside  the 
averaging  radius  is  sufficiently  rapid  to  produce  positive 
(sensible  only)  turbulent  heat  flux  from  the  atmosphere  to 
the  (Kean  as  SST  becomes  colder  than  air  temperature. 
This  impact  of  the  cyclones  is  likely  exaggerated  because 
the  air  temperature  used  to  force  the  model  is  fixed  and 
the  atmospheric  feedback  cannot  occur.  Some  atmo¬ 
spheric  feedback  is  present  in  the  COAMPS  air  temper¬ 
ature  forcing  fields  since  it  was  coupled  to  an  ocean  model, 
but  the  air  temperature  probably  remains  too  high  be¬ 
cause  the  coarse  representation  of  Ivan  in  the  CO  AM  PS 
atmospheric  model  presumably  generated  insufficient  SS'l 
cooling.  Experiment  GOM13  (no  H*W1ND  blending)  is 
consistent  with  this  hypothesis  since  it  underestimated  SST 
cooling  in  the  cold-core  cyclones  by  ^2°C  (Table  2).  De¬ 
spite  this  exaggerated  flux  reversal,  it  remains  clear  that 
thermal  energy  provided  by  the  ocean  to  the  inner-core 
region  of  Ivan  is  very'  sensitive  to  the  presence  of  these 
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FIG  9.  Air-sca  thermal  heat  Mux  averaged  over  (left)  I  9A\ia\ and  (right)  .vOA’max  Following  the  path  of  the  storm 
and  graphed  as  a  function  of  storm  central  latitude.  In  each  panel,  the  control  experiment  is  compared  to  alternate 
experiments  representing  a  particular  category  of  model  attribute  being  varied:  model  resolution  (top  row),  vertical 
mixing  choice  (second  row),  wind  stress  drag  coefficient  (third  row),  turbulent  heat  flux  drag  coefficients  (fourth  row ). 
and  both  surface  forcing  resolution  and  the  outer  model  used  for  nesting  (bottom  row).  The  control  experiment  is 
represented  by  the  black  line  in  all  panels.  Colors  representing  other  experiments  are  shown  in  the  legends  at  right 
Negative  values  indicate  heat  flux  from  the  ix*ean  to  the  atmosphere 


warm  and  cold  ocean  features.  The  progression  of  Q/over 
the  larger  3/?ma\  averaging  region  is  similar  over  the 
northwest  Caribbean  and  the  LC,  but  the  change  in  QT  i ts 
the  storm  passes  over  the  tv\o  cyclones  is  substantially 
reduced  because  warmer  water  outside  of  the  two  cyclones 
is  included  in  the  averaging  region.  The  flux  is  near  0  over 
the  southern  cyclone  and  remains  negative  (ocean  to  at¬ 
mosphere)  over  the  northern  one. 


Lillie  sensitivity  of  Q7  to  either  the  horizontal  or  ver¬ 
tical  grid  resolutions  exists  fur  both  averaging  radii  (Figs. 
9a  and  9b),  but  substantial  sensitivity  exists  for  all  other 
model  attributes.  Focusing  on  the  1  4 #max  averaging 
area,  altering  the  vertical  mixing  algorithm  produces 
differences  up  to  200  W  m  :  (20% -30%)  over  the  north¬ 
west  Caribbean  Sea  and  LC  in  the  GOM  (Fig.  9c).  Al¬ 
tering  either  the  wind  stress  or  turbulent  heal  flux  drag 
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Flo.  10.  Time  series  of  H^al  (Lfi )  a  location  in  l  he  eastern  part  of  the  southern  cold  cyclone  and  (right )  the  eastern 
side  of  the  warm  ring  at  the  two  locations  shown  in  the  lop  panel  of  I  ig.  5  I  lie  experiments  arc  sorted  into  separate 
panels  based  on  the  type  of  model  attribute  as  in  Fig.  9.  I  he  control  experiment  is  represented  by  the  black  line  in  .ill 
panels.  Colors  representing  other  experiments  arc  shown  in  the  legends  at  right. 


coefficients  produces  larger  differences  of  up  to  300  W  m 
(30% -40%)  for  CD  (Fig.  9c)  and  250  W  m  2  (25%- 
30%)  for  Cxi i  and  C\  s  (Fig.  9g).  Although  heal  flux 
sensitivity  to  changes  in  Cm  and  CYs  is  large,  SST  is  in¬ 
sensitive  since  surface  heat  flux  does  not  make  a  leading- 
order  contribution  to  the  SST  change  based  on  results 
from  earlier  studies  (c.g.,  Jacob  et  al.  2(XK);  Shay  and  Jacob 
2(XKi).  Large  sensitivities  arc  observed  as  expected  for  both 
the  low-resolution  forcing  and  the  alternate  model  ini¬ 
tialization  experiments  (Fig.  9i).  Except  for  turbulent  flux 
drag  coefficients,  changes  in  model  attributes  that  cool  SST 
faster  produce  reduced  heat  loss  to  the  atmosphere  as 
exacted 

Sensitivity  results  arc  generally  similar  over  the  larger 
3/?man  averaging  radius.  However,  over  the  northwest 
Caribbean  Sea  and  LC  in  the  GOM,  the  differences  that 
result  from  altering  the  vertical  mixing  and  wind  stress 
drag  coefficient  arc  smaller  than  the  differences  over  the 


smaller  averaging  radius  (Fig.  9d  and  9f).  The  sensitivity 
of  surface  heat  flux  to  altering  the  model  attributes  is 
greatest  in  the  inner-core  region  of  the  storm.  The  im¬ 
pacts  on  hurricane  intensity  will  therefore  depend  on  the 
poorly  known  radial  distance  within  which  the  ocean 
actually  provides  thermal  energy  to  the  storm  (Shell  et  al. 
2002;  Cione  and  Uhlhorn  2003). 


6.  Ocean  dynamics 

a.  Sensitivity  of  vertical  velocity  to  model  attributes 

Given  that  vertical  advcction  can  contribute  signifl 
cantly  to  SST  cooling,  the  sensitivity  of  vertical  advcction 
to  changes  in  he  model  attributes  is  assessed  by  graphing 
the  vertical  velocity  at  50-m  depth  (u's<>)  ns  a  function  of 
time  (Fig.  10)  at  the  same  two  locations  (warm  anticy¬ 
clone  and  southern  cyclone)  where  the  impacts  of  vertical 
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velocity  on  SST  cooling  were  analyzed  in  Fig.  6.  The  only 
difference  is  that  h'^o sampled  by  the  synthetic  instalments 
is  vertically  interpolated  to  50  m  instead  of  the  OML  base. 
The  behavior  of  vvso  m  the  control  experiment  is  consis¬ 
tent  with  the  results  summarized  in  Fig.  6.  with  maximum 
positive  n’5(,  occurring  near  the  time  of  Ivan's  passage  in 
the  southern  cyclone  and  with  somewhat  weaker  maxi¬ 
mum  positive  wso  lagging  passage  by  a  few  hours  along  the 
eastern  periphery  of  the  warm  anticyclone. 

The  interpretation  of  sensitivity  in  Fig.  10  is  confined  to 
the  lime  interval  when  direct  forcing  by  Ivan  was  im¬ 
portant  (roughly  ±3  IP).  The  least  sensitivity  in  u  soisduc 
to  the  altered  vertical  resolutions  in  both  ocean  features. 
Sensitivity  to  horizontal  resolution  is  small  except  in  the 
cyclone  where  )t  is  reduced  by  about  30%  at  the  lime 
of  maximum  forcing.  Little  sensitivity  is  evident  when 
Cel  and  Cp s  are  changed  as  expected  due  to  the  lack  of 
SST  sensitivity.  When  other  model  attributes  are  altered, 
maximum  ivs<>  can  vary'  by  30%  or  more  in  the  southern 
cyclone.  This  sensitiv  ity  is  sufficiently  large  that  differences 
in  the  upwelling  response  beneath  the  inner  core  of  Ivan 
could  contribute  significantly  to  differences  in  SST  cooling. 
By  comparison,  substantially  less  sensitivity  is  observed  in 
the  warm  ring  for  all  model  attributes.  In  this  ocean 
feature,  a  sharp  peak  in  upwelling  velocity  is  observed  0. 1 
IP  after  storm  passage  for  the  two  experiments  that  used 
the  uncapped  Large  and  Pond  and  the  Shay  and  Jacob  Cn 
parameteri/ations,  presumably  related  to  the  fact  that 
these  choices  have  large  values  at  high  wind  speed.  A 
secondary  peak  is  observed  at  0.3  IP  for  the  Shay  and 
Jacob  choice,  which  is  possibly  related  to  the  fact  that  this 
choice  has  comparatively  large  values  at  the  intermediate 
wind  speeds  found  in  the  outer  reaches  of  the  storm. 

b.  Velocity  profiles  produced  by  the  control 
experiment 

ADCP  observations  from  SEED  moorings  1, 3,  7,  and  9 
are  used  to  evaluate  the  capability  of  the  control  experi¬ 
ment  to  simulate  the  wind-driven  near-incrtial  current 
response  (Teague  et  al.  2007).  which  is  critically  important 
to  the  shear-driven  turbulent  entrainment  at  the  OML 
base.  The  u  and  u  components  at  these  moorings  are 
shown  for  the  upper  150  m  of  the  water  column  over  a 
time  interval  ranging  from  0.5  IP  before  the  storm  to  1.5  IP 
after  passage,  which  represents  the  forced  response  and 
the  subsequent  transition  to  near  inertial  current  variabil 
ilv  (Fig.  11).  Observed  velocity  components  are  graphed 
over  the  same  time-depth  range  except  for  the  near-surface 
region  where  the  ADCP  could  not  sample  velocity  because 
ofsidclobc  inlcrfcrenee.  The  responses  at  moorings  7  and  9 
where  the  water  depths  (>500  m)  arc  much  greater  than 
the  OML  thickness  are  dominated  by  forced,  near-incrtial 
waves,  and  the  simulated  response  closely  resembles  the 


observed  response  over  their  common  depth  range.  The 
magnitude  of  the  response  at  mooring  7  is  reduced  by 
the  unfavorable  anticlockwise  rotation  of  the  wind  stress 
vector  as  Ivan  passes.  ITiese  observations  and  simulations 
suggest  that  vertical  energy  propagates  out  of  the  surface 
mixed  layer  and  into  the  thennocline  in  a  manner  consis¬ 
tent  with  theory.  At  both  mtxirings,  westward  mean  flow  is 
present  in  the  //  held  while  the  mean  v  is  near  0. 

Although  the  near  inertial  time  scales  of  the  response  al 
shallow- water  (50  m)  moorings  1  and  3  are  generally  ev¬ 
ident,  the  response  is  more  complicated  than  in  deeper 
water.  The  v  response  al  mooring  1  displays  strong  mid 
depth  northward  How  between  0.5  and  1  IP  following 
storm  passage  that  is  captured  in  the  model  current  re 
sponsc.  The  observed  u  and  v  responses  at  mooring  3  re 
veal  a  relatively  strong  northeastward  flow'  event  in  the 
bottom  10  m  around  0.5  IP  after  storm  passage  that  is  not 
captured  by  the  model.  Without  observations  of  the 
three-dimensional  ocean  structure  in  the  vicinity  of  these 
moonngs.  we  cannot  determine  the  reason  for  this  near- 
bottom  How  event  or  why  the  model  did  not  reproduce  it. 
Overall,  however,  the  resemblance  between  the  obser¬ 
vations  and  simulations  al  these  moorings  is  encouraging, 
particularly  for  the  deeper  moonngs 

r.  Sensitivity  of  velocity  profiles  to  mode I  attributes 

A  Taylor  diagram  is  used  to  quantify  the  sensitivity  at 
SEED  mooring  9  by  comparing  the  v  fields  from  each  al¬ 
ternate  experiment  to  the  v  generated  by  the  control  ex 
periment  (Fig.  12.  top).  To  const  met  this  diagram,  the  mean 
was  first  removed  from  each  r  field  over  the  time  interval 
between  0.5  and  1.5  IP  and  over  the  upper  150  m  The 
resulting  two-dimensional  fields  were  normalized  by  the 
variance  of  the  held  from  the  control  experiment.  In 
contrast  to  the  situation  farther  offshore,  the  alternate 
initialization  (GOMI4)  is  noi  an  ouilicr  for  simulated 
velocity.  The  initial  flow  at  SEED  mooring  9  is  westward 
in  both  GOM14  and  GOMl  (Fig.  1 1),  so  the  different 
initializations  produced  similar  background  flows  in  this 
region. 

The  greatest  outliers  in  the  mooring  9 sensitivity  analysis 
are  the  low-resolution  atmospheric  forcing  experiment 
GOM13  and  the  one-dimensional  experiment  GOM15.  In 
GOMl 3,  velocity  amplitude  is  underestimated  by  up  to 
40%  relative  to  the  control  experiment.  TTie  correlation 
with  this  control  experiment  is  also  reduced  to  0.9,  pn 
marily  because  the  weaker  forced  currents  do  not  penetrate 
as  deep  into  the  water  column  (not  shown).  In  GOM15, 
velocity  amplitude  is  overestimated  by  >40%  and  the 
correlation  with  the  control  experiment  is  reduced  to  0.66 
because  the  response  is  purely  inertial  and  energy  does  not 
penetrate  into  the  ocean  interior  through  near  inertial 
wave  dynamics  (Shay  and  ELsherry  1987:  Shay  et  al.  1989) 
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The  correlation  is  also  reduced  because  GOM15  produces 
current  fluctuations  with  an  exactly  inertial  period  while 
GOM1  produces  fluctuations  with  a  period  about  10% 
smaller  than  inertial  (Fig.  1 1 ).  The  impacts  of  upweJIing  on 
OML  thickening  arc  also  absent  in  GOMI5, 

The  next  greatest  level  of  sensitivity  results  from  alter¬ 
ing  the  wind  stress  drag  coefficient.  The  choices  that  pro¬ 
duce  larger  (smaller)  Cn  values  produce  larger  (smaller) 
RMS  amplitudes  than  the  control  experiment.  Both  in¬ 
creased  and  decreased  Cp  are  associated  with  a  reduction 
in  correlation,  primarily  because  increased  (decreased)  C/> 
is  associated  with  increased  (decreased)  penetration  of  the 
wind-forced  currents  deeper  into  the  water  column  (not 
shown).  Two  choices  of  C»  (Large  and  Pond  capped, 
GQM9;  Jarosz  el  al.,  GOM1 1 )  display  little  difference 
from  the  control  experiment.  These  results  are  gener¬ 
ally  consistent  with  simulations  of  the  ocean  response 
to  Hurricane  Francis  (2004)  by  Sanford  et  al.  (2007), 
who  determined  that  runs  with  CL>  values  capped  or  re¬ 
duced  at  high  wind  speeds  outperformed  runs  using  the 
uncapped  Large  and  Pond  representation.  Low  sensitiv  it\ 
is  observed  with  altering  vertical  mixing  choice  and  with 
both  vertical  and  horizontal  resolutions.  The  alternate 
heat  flux  coefficient  experiment  produces  fluctuations  that 
arc  nearly  identical  to  the  control  experiment,  as  expected. 

d.  Evaluation  of  velocity  profiles  against 
observations 

A  Taylor  diagram  analysis  is  performed  to  compare  v 
fluctuations  at  SEED  mooring  9  produced  by  all  15  ex¬ 
periments  to  the  observed  fluctuations  (Fig.  12.  bottom) 
It  is  conducted  over  the  common  depth  range  of  40- 
150  m  and  from  0.51  P  before  to  1 .51 P  after  the  storm.  All 
fields  are  normalized  by  the  variance  of  the  observed  v 
fluctuations.  This  analysis  is  supplemented  bv  skill  scores 
calculated  over  the  same  depth-time  range  (Table  3).  The 
control  experiment  produces  one  of  the  most  realistic 
simulations  of  upper-ocean  v  fluctuations  with  (eorrela 
lion,  normalized  RMS  amplitude,  normalized  RMS  dif 
fere  nee)  of  (0.89,  1.46.  0.48)  and  with  skill  scores  of  0.71 
and  0.77  for  u  and  v%  respectively. 

Vertical  resolution,  the  turbulent  heat  flux  drag  co¬ 
efficient.  and  both  the  capped  Large  and  Pond  and  Jarosz 
cl  al.  choices  of  wind  stress  drag  coefficient  all  have  only 
a  minor  impact  on  the  quality  of  the  response  while  hor¬ 
izontal  resolution  has  only  a  slightly  larger  influence.  ’T  he 
vertical  mixing  choice  has  a  substantial  influence  on  the 
magnitude  of  the  current,  with  both  the  MY  and  GISS 
experiments  (GOM5  and  GOM6)  producing  a  substan¬ 
tially  weaker  velocity  response.  This  result  contrasts  with 
the  sensitivity  evaluation  in  Fig.  12  (top)  where  the  other 
vertical  mixing  choices  produced  velocity  fluctuations  with 
magnitudes  close  to  that  of  the  control  experiment.  This 


difference  occurred  because  the  present  comparison  to  the 
observations  is  limited  to  the  depth  range  between  40  and 
150  m.  suggesting  that  the  forced  velocity  fluctuations  in 
GOM5  and  GOM6  are  more  surface  intensified  and  do 
not  penetrate  downward  into  the  water  column  as  effec¬ 
tively  as  they  do  in  the  control  experiment.  Consequently, 
the  KPP  mixing  choice  produces  a  more  realistic  current 
response  with  respect  to  these  observations.  The  wind 
stress  drag  coefficient  has  a  significant  impact  on  the  pre 
dieted  velocity  profile.  The  choices  of  C/>  other  than  I xirge 
and  Pond  capped  and  Jarosz  et  al.  degraded  the  overall 
realism  of  the  simulated  velocity  fluctuations.  Tlie  two 
para  me  l  c  riza  lions  that  continue  to  increase  in  value  al  high 
wind  speeds  produced  less- realistic  velocity  responses,  in 
agreement  with  Sanford  et  al.  (2007).  The  two  most  un 
realistic  simulations  (GOM13  and  GOM15)  were  poorly 
correlated  with  the  observations  and  had  low  skill  scores 
(Table  3). 

Taylor  diagrams  comparing  the  simulated  v  responses 
to  the  observations  are  also  presented  for  moorings  1  and 
7  (Fig.  13).  In  general,  the  conclusions  reached  from  u  at 
mooring  9  (Fig.  12,  bottom)  are  upheld,  in  particular  the 
dependence  on  the  choice.  Comparisons  arc  again 
poor  for  low -resolution  forcing  (GOM13)  and  for  one 
dimensional  ocean  dynamics  (COM  15).  For  the  most 
realistic  experiments,  the  qualitv  is  slightly  reduced  in 
comparison  to  mooring  9  due  to  the  weaker  response  to 
the  west  of  the  eye.  A  further  reduction  in  quality  is  ev  ¬ 
ident  over  the  shelf  with  maximum  correlations  close  to 
0.7  as  opposed  to  between  0.8  and  0.9  farther  offshore. 

7.  Conclusions 

The  response  of  an  ocean  model  to  Hurricane  Ivan 
(2004)  has  been  documented  in  the  GOM  and  NW  Ca¬ 
ribbean  Sea.  The  sensitivity  ol  the  simulated  response  to 
changes  in  several  specific  attributes  of  the  model  (Table 
I )  was  quantified  and  the  realism  of  the  model  response 
was  evaluated  against  SSI  images  derived  from  satellite 
and  in  situ  measurements  along  with  upper-ocean  current 
profiles  measured  at  ADCP  moorings.  Key  results  are 
summarized  in  Table  4,  along  with  recommendations  for 
each  model  attribute.  In  drawing  conclusions  from  these 
results,  it  must  be  kept  in  mind  that  they  are  represent 
tive  of  one  particular  storm  in  one  distinct  ocean  region, 
and  have  been  obtained  from  a  stand  alone  ocean  model 
without  atmospheric  feedback  that  was  driven  by  quasi¬ 
optima  1  but  imperfect  atmospheric  forcing.  Despite  these 
limitations,  tlie  present  study  represents  a  baseline  effort 
to  understand  the  strengths  and  limitations  of  ocean 
models  with  respect  to  their  numerics  and  paranieleriza- 
lions.  a  first  step  toward  the  ultimate  goal  of  improving 
their  performance  in  coupled  forecast  models  ll  provides 
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guidance  toward  where  we  should  most  effectively  tnvest 
our  efforts  to  improve  model  performance. 

'Hie  present  results  demonstrate  that  representation 
of  the  ocean  by  one-dimensional  ocean  models  may  be 
inadequate  for  predicting  the  ocean  response  to  some 
storms.  In  the  region  studied  in  this  analysis,  energetic 
ocean  features  sueh  as  the  LC  and  eddies  distort  the  wind- 
driven  momentum  response  and  the  associated  upwelling 
patterns  that  can  significantly  influence  SST  cooling  be¬ 
neath  a  storm  (e.g.,  Jaimes  and  Shay  2009,  2010).  This 
impact  was  particularly  evident  between  the  detaching 
warm  anticyclone  and  a  cold-corc  cyclone  present  just 
to  the  southeast  of  this  ring  where  upwelling  was  en¬ 
hanced  within  the  cyclone.  However,  limitations  of 
one-dimensional  ocean  models  may  still  be  significant  in 
less-energetic  ocean  regions,  particularly  for  relatively 
slow-moving  storms  in  regions  with  a  relatively  thin  sur 
face  warm  layer.  The  absence  of  upwelling  may  signifi¬ 
cantly  reduce  the  SST  cooling  in  these  cases,  f  urthermore, 
one-dimensional  ocean  models  predict  a  purely  inertial 
current  response  with  energy  that  cannot  propagate 
into  the  ocean  interior  through  wave  dynamics  These 
factors  may  be  especially  important  when  a  storm  is 
closely  followed  by  a  second  one  (c.g.,  Isidore  and  Till 
in  2002;  Katrina  and  Rita  in  2005;  Gustav  and  Ike  in 
200S). 

The  urge  to  design  efficient  coupled  forecast  models 
containing  a  one-dimensional  ocean  model  is  understand¬ 
able  given  operational  eonstraints  from  available  computer 
resources.  Although  the  choice  of  vertical  mixing  pa¬ 
rameterization  has  little  effect  on  run  time,  it  was  reduced 
by  about  70%  in  the  one-dimensional  experiment.  Given 
the  results  of  this  study,  how  ever,  we  cannot  conclude  that 
a  one-dimensional  model  will  always  be  adequate  for 
performing  intensity  forecasts.  Furthermore,  the  western 
boundary  regions  of  ocean  basins  where  the  use  of  three- 
dimensional  ocean  models  is  most  likely  to  be  important 
are  located  close  to  land  where  intensity  forecast  er¬ 
rors,  particularly  in  regard  to  rapid  intensification  and 
rapid  weakening,  will  have  the  worst  consequences.  These 
questions  will  ultimately  have  to  be  answered  by  coupled 
forecast  experiments  for  a  large  number  of  storms  in  dif¬ 
ferent  ocean  regions  and  basins. 


Existing  operational  atmospheric  models  are  run  at  a 
horizontal  resolution  that  is  loo  coarse  to  resolve  the  inner 
core  TC  structure,  overestimating  WMAX  by  a  factor  of 
2  and  underestimating  the  maximum  wind  by  40%- 
50%  in  the  case  of  the  27-km  COAMPS  model  For  this 
reason,  the  H*W1ND  wind  analysis  was  blended  into 
fields  from  the  U.S.  Navy  CO  AMPS  atmospheric  model 
to  generate  forcing  for  the  control  experiment  in  this 
study  that  at  least  marginally  resolved  the  horizontal 
scales  of  the  eye  and  eycwall.  An  alternate  experiment 
forced  by  the  CO  AM  PS  model  alone  substantially  de¬ 
graded  the  quality  of  the  simulated  ocean  response.  Re¬ 
search  studies  of  the  ocean  response  to  a  TC  should  not 
use  forcing  fields  from  these  existing  operational  atmo 
spheric  models.  The  question  of  whether  the  scales  re 
solved  by  H*W1ND  are  fully  adequate  could  not  be 
addressed  here  because  higher-resolution  forcing  from  an 
atmospheric  model  that  accurately  resolved  both  the 
track  and  intensity  of  Ivan  was  not  available.  Since  op¬ 
erational  forecast  centers  arc  moving  toward  using  nested 
atmospheric  models  with  resolution  of  0(1  km),  atmo¬ 
spheric  model  resolution  should  not  be  an  issue  for  cou 
pled  TC  forecast  models  in  the  future. 

The  accuracy  of  the  upper  ocean  thermal  response  to 
Ivan  was  most  sensitive  to  the  ocean  model  initialization. 
The  two  important  aspects  of  initialization  are  ocean 
feature  location  and  providing  accurate  profiles  of  tem¬ 
perature,  salinity,  and  densit\  within  the  existing  features. 
The  control  experiment  was  initialized  from  a  data 
assimilative  ocean  hindcast  where  I  he  LC.  a  detaching 
warm  ring,  and  two  associated  cold-core  cyclones  that  had 
a  large  impact  on  the  final  SST  cooling  pattern  were 
correctly  situated  An  alternate  experiment  initialized 
from  a  non-assimilative  ocean  model  demonstrated  that 
feature  location  had  a  much  larger  influence  on  the  re 
sponses  of  SST,  wind-driven  upwelling.  and  heat  flux  front 
the  ocean  to  the  atmosphere  than  did  other  model  and 
air  -sea  flux  parameterizations.  Although  this  study  w  as 
conducted  in  a  region  where  water  mass  and  heat  content 
differences  between  ocean  features  was  large,  large  biases 
in  initial  temperature,  salinity,  and  density  profiles  along 
with  the  thickness  of  the  warm  layer  can  also  significantly 
degrade  upper-ocean  SST  forecasts  in  all  ocean  regions 


Fig  12.  (top)  T«ylor  diagram  comparing  the  simulated  north  -south  (cross  shelf  and  along  track)  velocity  component  from  lhe  control 
experiment  to  the  same  field  from  all  of  lhe  remaining  experiments  over  the  upper  150  m  ai  SEED  mooring 9  Symbols  and  colors  follow 
the  conventions  used  in  the  lop  panel  of  Fig.  <S.  The  large  black  square  represents  a  perfect  comparison.  (|vot tom)  Taylor  diagram 
comparing  lhe  simulaled  north  south  (cross  shelf  and  along  l rack )  velocity  component  for  all  experiments  lo  the  same  field  from  the 
ADC'P  observalions  at  SEED  moormg  9.  I  he  comparison  was  performed  over  the  common  depth  range  where  hoili  observed  and 
simulaled  data  were  available  within  the  upper  I  SO  m  of  lhe  water  column,  as  shown  in  Fig  1 1  Symbols  and  colors  follow  the  conventions 
used  in  the  hotiom  panel  of  Fig  S'.  The  large  black  square  again  represents  a  perfect  comparison 
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I’ABt  I.  4.  Summary  of  lhc  impacts  of  the  changes  in  different  model  attributes  on  the  SS  I  response  pattern,  the  heal  flux  from  the  ocean 
to  the  atmosphere,  the  wind-driven  upwelling,  and  the  upper  ocean  near-inerti.il  momentum  response  Heal  flux  differences  are  for  the 
flux  averaged  over  l  4  /?MAX.  Sensitivity  is  measured  relative  to  the  control  experiment  If  the  sensitivity  depends  on  the  location  relative  to 
the  energetic  ocean  features,  the  impacts  on  sensitivity  are  modified  by  the  word  “potentially”  to  indicate  that  this  level  of  sensitivity  may 
not  be  true  everywhere. 
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on  the  turbulent  heal 

high  impact  on  all  important  aspects  of  model  sensitivity  flux  as  expected,  but  had  very  liltle  impact  on  SST  cooling 

(Table  4).  For  wind  stress,  the  magnitude  of  Cn 

al  high  since  entrainment  at  the  OMI  bast 

makes  the  dominant 

wind  speeds  impacts  SST  cooling  and  heat  flux  from  the  contribution  to  this  etxilitig  (e.g.,  Jacob  et  al.  2000;  Shay 

ocean  to  the  atmosphere  by  1 )  modifying  the  rate  of  OMI  and  Jacob  2006).  Substantial  effort  w 

ill  hav e  to  be  directed 

deepening  by 

wind-driven  turbulence  (friction  velocity)  toward  verifying  the  parameterizations  of  the  surface  flux 

and  2)  altering  the  magnitude 

and  three-dimensional  drag  coefficients  in  coupled  TC  forecast  models. 

structure  of  wind-driven  currents  and  wind-driven  upwell-  Although  the  SST  response  displayed  intermediate 

tng.  For  turbulent  heat  flux,  an  alternate  experiment  using  sensitivity  to  the  three  vertical  mixing  models  tested,  both 

I  Hi  )3  Taylor  diagrams  comparing  the  simulated  north  south  (cross  shelf  and  along  track  )  velocity  component  for  all  experiments 
to  the  same  field  from  the  ADC  P  observations  at  SFFD  moorings  (lop)  1  and  (bottom)  9.  The  comparison  was  performed  over  the 
common  depth  range  where  both  observed  and  simulated  data  were  available  within  the  upper  I  s()  m  nt  the  vv alcr  column,  as  shown  in 
Fig.  1 1  Sy liiHols  and  colors  follow  the  conventions  used  in  the  bottom  panel  of  f  ig.  N  The  large  black  square  represents  a  perfect 
comparison. 
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wind-driven  up  welling  and  air-sea  heat  flux  displayed  high 
sensitivity  that  was  comparable  to  the  high  sensitivity  to 
surface  flux  drag  coefficients.  Although  a  substantial  effort 
to  improve  the  performance  of  vertical  mixing  models  is 
warranted,  further  improvement  in  model  initialization 
and  surface  flux  parametcrizations  will  first  be  necessary': 
otherwise,  vertical  mixing  may  be  “tuned"  to  correct  for 
errors  and  biases  resulting  from  these  other  factors. 

Horizontal  resolution  has  a  low  to  intermediate  in 
fluence  on  all  important  aspects  of  model  performance 
(Table  4),  Although  it  has  intermediate  impact  on  VSST 
cooling,  it  has  only  a  small  influence  on  the  heat  flux  from 
the  ocean  to  the  atmosphere.  This  probably  results  in  part 
because  much  of  the  RMS  difference  (X'curs  at  small  scales 
since  the  higher-resolution  experiment  provides  a  sharper 
resolution  of  temperature  changes  across  fronts.  As  a  re¬ 
sult  these  SST  differences  have  only  a  small  impact  on 
area  integrated  heat  flux.  Sensitivity  to  vertical  resolution 
is  low  with  respect  to  all  important  aspects  of  model  per¬ 
formance.  Given  the  greater  sensitivity  to  other  model 
attributes,  there  is  no  need  to  invest  a  lot  of  effort  into 
optimizing  the  model  s  horizontal  or  vertical  resolution 
Based  on  the  Ivan  results,  ocean  models  run  at  a  horizontal 
resolution  of  <10  km  and  a  vertical  resolution  within 
the  OML  of  <10  m  should  be  adequate  for  addressing 
the  greater  concerns  enumerated  above  since  further 
increases  in  resolution  will  have  a  small  impact  on  SST 
cooling  and  on  the  heat  (lux  from  the  ocean  to  the 
atmosphere. 

Although  these  results  are  based  on  one  ocean  model 
for  one  hurricane  in  one  ocean  region,  they  represent  an 
early  step  toward  designing  strategies  that  will  improve 
the  performance  of  the  ocean  component  of  coupled 
forecast  models.  Ocean  model  studies  must  be  extended 
to  other  TCs  and  other  ocean  regions,  while  evaluation 
and  sensitivity  analysis  must  also  be  extended  to  coupled 
forecast  models.  These  efforts  will  require  ongoing  pro¬ 
grams  to  obtain  high  quality  observations  of  currents  and 
baroclinie  structure  with  sufficient  detail  and  both  three- 
dimensional  and  temporal  patterns  of  coverage  to  per¬ 
form  the  detailed  evaluation  sludies  requtred  to  formulate 
effective  model  unprovcmcnl  strategies. 
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